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Summary. We consider the continuum version of an elastic chain supporting 
topological and non-topological degrees of freedom; this generalizes a model for 
the dynamics of DNA recently proposed and investigated by ourselves. In a 
certain limit, the non-topological degrees of freedom are frozen, and the model 
reduces to the sine-Gordon equations and thus supports well-known topological 
soliton solutions. We consider a (singular) perturbative expansion around this 
limit and study in particular how the non-topological field assume the role of 
an auxiliary field. This provides a more general framework for the slaving of 
this degree of freedom on the topological one, already observed elsewhere in 
the context of the mentioned DNA model; in this framework one expects such 
phenomenon to arise in a quite large class of field-theoretical models. 

1 Introduction 

In a recent paper [1] we analyzed a "composite" model of DNA torsion dynamics, 
pretty much in the spirit of the by now classical Peyrard-Bishop, Yakushevich 
and Barbi-Cocco-Peyrard models [2, 3, 4, 5, 6, 7, 8, 9, 10]; this model amounted 
to a field-theoretic Lagrangian with topological and non-topological interacting 
fields. We studied in particular the solitary wave excitations (we call these 
solitons for ease of language) it supports; this analysis displayed some puzzling 
and somehow surprising features: 
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(a) on the one hand, the speed of given sohton solutions of a certain (relevant) 
type is not a free parameter but turns out to be fixed by (the parameters 
of) the model; 

(h) On the other hand, when performing a perturbative analysis (exact solu- 
tion of the model was not possible) the non-topological field turns out to 
be completely determined - we say then it is slaved - by the topological 
one. 

It appeared that these phenomena are rather obviously not specific to 
the model considered there, nor to DNA dynamics, but actually common to 
a much wider class of models; they coiild also be of rather obvious interest 
in applications. We believe therefore they arc worth further investigation and 
clarification. 

In a previous related paper [11] we have investigated the selection of the 
soliton speed, i.e. point (a) above; in the present paper we focus on point (b), 
i.e. on the slaving of one field. We identify the origin of this in the remarkable 
fact that the slaved field, while being on equal footing in the full theory, becomes 
an auxiliary field in the perturbative expansion; that is, at each order in the 
perturbative expansion the Lagrangian depends on the field but not on the 
conjugate momentum.^ 

Auxiliary fields appear in a number of field theories (e.g., auxiliary fields 
play a crucial role in supersymmetric field theories where they allow for off-shell 
closure of the algebra [13, 14]), but it is somehow surprising to have these in the 
framework of elastic chains, and even more so in the context of DNA dynamics. 
Moreover, the appearance of auxiliary fields and the; rc;latc;cl slaving mechanism 
are quite interesting from the point of view of the would-be functional role of 
DNA solitons, in that they would provide a robust mechanism for the coordina- 
tion of different degrees of freedom, essential if the solitons have to play a role - 
as conjectured since a long time [15] (see also [4, 5, 9]) - in DNA transcription. 

More generally, the mechanism at work here can be described as follows: 
one considers a system with two degrees of freedom, and the (singular) limit in 
which one of the two fields is constrained to zero, determining an exact solution 
to the full field equations for the other field; if one considers the (singular) 
perturbation expansion around this field configuration, it turns out that the 
perturbation for one field are completely determined - via algebraic identities 
- by the solutions for the dynamical equations for the other. It appears the 
mechanism is in these abstract terms - quite general; we expect therefore that 
it can be generalized (as the speed selection mechanism mentioned above [11]) 
to a number of other physical situations. 

In order to show more clearly the mechanism at work, we study an ele- 
mentary model, made of a chain of coupled double pendulums; while the first 

'^From the point of view of DNA dynamics it is interesting to remark, in tiiis respect, that 
the BCP model of DNA dynamics was recently shown [12] to admit a constant of motion in its 
field-theoretic limit; the presence of a constant of motion leads to a reduction in the effective 
dimensionality of the theory, similarly to what happens in the presence of an auxiliary field. 
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Figure 1: Notation for each of the double pendulums along the chain; see text 
for detail. 

pendulums are standard ones, free to swingle round the circle, the second pendu- 
lums oscillations are constrained to an angular range \ip\ < (po ■ We consider the 
continuum limit of the chain, which yields a field-theoretic Lagrangian with a 
topological and a non-topological degree of freedom (corresponding respectively 
to the angular coordinates of first and second pendulums). 

2 The double pendulums chain model 

We will consider an infinite chain of double pendulums, suspended at points 

= n6 on a straight horizontal line, denoted below as the "axis" of the chain. 

Each of the double pendulums is made of a first massless rigid beam of 
length R suspended at one extremum on the line and having a point mass M at 
the other extremum; and of a second massless beam of length r, joined at one 
extremum with the extremum of the first beam where the point mass is, having 
a point mass m at the other extremum. 

The pendulums can only rotate in the plane orthogonal to the axis. The 
rotation angle of the first pendulum at site i will be denoted as Of, the angle 
of rotation of the second pendulum (with respect to the direction of the first 
pendulum) will be denoted as Sec fig.l. 

The second pendulum is not free to swingle through a full circle, but is 
instead constrained to stay in the range < 4>q. This constraint will also be 
modelled by adding a constraining potential, i.e. a potential Vc(</>) which has the 
effect of limiting de facto the excursion of the (j) angles. The rest configuration 
will correspond to Oi = (pi = Q. 

Let us denote, dropping for a moment the index i, as {X^Y) and {x,y) the 
cartesian coordinates - in the plane orthogonal to the double helix axis ^, with 
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Figure 2: Pendulums at successive sites k and k + 1 interact via harmonic 
potentials coupling the masses at corresponding positions on the two pendulums; 
see text. 



origin in the suspension point - of the point mass at the end of the first and 
of the second pendulum respectively. The vertical direction will be along the x 
axis. 

In the equilibrium position, given hy 9 — (j) — 0, we have (X, Y) = {R, 0) 
and {x, y) = {R + r, 0). In general, it will be 

X = Rcos6 , X = Rcos9 + rcos{9 + (I)) ; , . 

Y = Rsm9 , y = Rsm9 + rsm{9 + (l)) . ^' 

The kinetieal energy for the double pendulum at site i will then be Tj = 
{l/2)[M{{Xif + {Yi)^)+m{{xi)^ + {yif)]-, using (1) to express , and Xi,yi 
in terms of the angles 9i and we get 



J-i — 2 ' 2 



RHf + 2Rrcos^i (^{9if + 9i^,^ + r^{9i + ^if'\ . (2) 



The total kinetic energy is of course T = ^ . Tj. 

As for interactions, each pendulum will be coupled to nearest neighbors via 
harmonic potentials, i.e. linear springs attached to the M and to the m masses 
of the pendulums at sites k and fc + 1. Moreover, the mass m at the end of the 
second pendulums will experience a (site-independent) external potential. 

Thus we have the following forces, and correspondingly potential energy 
terms in the Lagrangian: 

(1) A coupling between successive "first pendulums", described by a potential 

= Kt[l - cos{9,+i - 9i)] ; (3) 



this will be referred to as the "torsional interaction" . The total torsional poten- 

j(a,i) 



tial is of course Vt = ^ [/^ 



(2) A coupling between successive "second pendulums" on each chain, described 
by a potential 

this will be referred to as "stacking interaction". In terms of the 6i and <j)i 
angles, using again (1) and with A = r/i?, we get 

C/i'^ = {ksR^/2) [{{cosei - cos6i+i) + A[cos(<?ii + Oi) - cos{4>i+i + (9i+i)])^ 
+ ((sin^i - sin^i+i) + X[sm{(f)i + 6i) - sin(^i+i + ^i+i)])^ ; 

(4) 

the total stacking potential is of course Vs = ^ v}"''^^ 

(3) Interaction with the external potential, 

C/« = F{Oi,^i) . (5) 

The total external potential will of course be = J2i Up^ ■ In our case, the 
external potential will just be the gravitational one, i.e. 

F{ei,^i) = g[M R{1- cos Oi) + m[R + r-{R cos ei +r cos{(j)i + ei))]] . (6) 

(4) Interaction with the confining potential, 

[/« = H<Pi) , (7) 

with h{x) a convex even function, essentially flat for <C ^O) and rising sharply 
at I a; I « 0o- The total external potential will of course 

We have thus provided explicitly all terms appearing in the total Lagrangian 

L = T - {Ut + Us + Up) . (8) 

The dynamics of the model will be governed by the corresponding Euler- 
Lagrange equations, 

dL d [dL\ 

" dt j = ' ' 

We will not write these explicitly; they are rather involved, due mainly to the 
kinetic term.^ 



^It would be possible to obtain a simpler form for most terms by using different coordinates, 
i.e. passing to {6, ip) with the angle made with the rest (vertical downward) direction by the 
line from the suspension point of the first pendulum to the mass of the second pendulum. This 
would however introduce a very complex dependence of the confining potential h{cp) = h{0, tp), 
which we prefer to avoid. 
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3 Continuum version 



If one is mainly interested in solutions varying little on the intcrsite scale, i.e. 
such that at all times — Oi\ and — are small - say of order e -C 1, 
we can pass to a continuum description. 

In this, the arrays 9i{t) and (l)i{t) will be replaced by fields 6(a;, t) and ^{x, t) 
such that, say, <^{kS,t) = (t>k{t) and e{kS,t) = Ok{t). 

Expanding in a Taylor series up to order two in 6 we have 



^x ±6,t) = t) ± 6<^^{x, t) + (5V2)$^^(a;, t) , 

Six ±s,t) = e(x, t) ± se^^ix, t) + {s^/2)e,^{x, t) . 



2/o^c> .^ (10) 



Inserting this into the Euler-Lagrange equations (9), we get the field equa- 
tions, see (12) below, governing our system in the continuum approximation. In 
order to write these, it is convenient to set 

a($) := 1 + (ii/r) cos $ , /3($) := (1 + (i?/r)^ + 2(E/r) cos $) (11) 

we will moreover set 

Kg := Kg (5^ , Kt := KfS'^ . 
With these, we get the field equations 

e,4Kt + Ksr^f3m + $.,[i^.r2a($)] = 

-rR[mM^t + 26*) - K^^^i^^ + 26^)] sin$+ 

+g[Rim. + M) sin 6 + mr sin($ + 6)] ; 

+/i'($) + rR{me1 - KsQl) sin $ + gmr sin($ + 6) . 



(12) 



We will from now on work in the continuum approximation, i.e. with (12). 
These are also obtained by rewriting the Lagrangian L in the continuum 
approximation; the explicit expression of this (written for future reference) is 

C = (1/2) [MR^Ql - KtQl + m {r\^t + Qt? + R^^l 
+2ri?et($t + et)cos$) 

+2g[{m + M)Rcose + mrcos{'t + Q)] -2h[<^] 
-K,[{RQ^ cos e + r{^x + ^x) cos($ + 9))' 
+ (i?e^sine + r($^ + e^)sin($ + e))'] . 

Then (12) are also obtained as the associate Euler-Lagrange equations 
„ dL d f dL\ d f dL\ „ 

^* d^-di[wj-d^[M:)=^- ^^^^ 
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The field equations (12) should be supplemented by a specification of the 
function space the solution are required to belong to; these can be given as 
boundary conditions. 

The physically natural condition is that of finite energy at any time t; that 
is, considering the energy density 

H{x,t) = H[e{x,t),^x,t)] = T+{Ut + U, + Up + U,) , (15) 

and fixing the additive arbitrary constant so that the minimum of the potential 
energy is at zero, the requirement that 

/ + 00 
H{x,t)dx < 00. (16) 
-oo 

For this to be satisfied, in view of the form of T, we must require that 

lim = lim Ba; = , (17) 

x— >iboo X— >±oo 

i.e. that <i> and O are asymptotically constant in x. 

Moreover, the fields $ and should go at minima of the potential energy for 
X — » ±oo; in view of the form of the potential energy, we must actually require 
that 

lim $ = ; lim e = 27rn± . (18) 

(Note that the difference between topological and non topological degrees of 
freedom shows up here.) 

We can always change origin of the 6 angles, so that 

lim e = , lim e = 27riV . (19) 

X— * — oo X— > + oo 

Thus, finite energy solutions possess a topological index, the integer N, 
identified by the asymptotic behavior at x ±oo. One gets easily convinced 
that these boundary conditions (and the finite energy condition) are preserved 
under the time evolution described by (12). 

4 Travelling wave solutions 

We are specially interested in travelling wave (TW) solutions, i.e. solutions such 
that 

&{x, t) = ^{x - vt) = 'diz), $(x, t) = ip{x - vt) = ip{z) . (20) 

It is immediate to see that the boundary conditions implied by the require- 
ment of finite energy become in this framework 

lim^^ioo f'iz) = , lim^^ioo i^'iz) = ; .^i) 
lim^^ioo ^{z) = , lim^^ioo "^iz) = 27rn± . 

We will again choose variables so that ^?(— oo) = 0, i?(+oo) = 27rN. 
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Inserting the ansatz (20) into the equations (12), and writing M = pm for 
ease of notation, we get the equations for traveUing wave solutions: 

fJ-^zz {r^ + rR cos if) + [Kt - K.R^p + ii{r^ + R^ + R^p + 2rR cos <^)] + 
— /xr-R^^ sint^ — 2/iri?<^2i?2 sint^ — rn5[i?(l + p) sini? + r sin((/7 + 1?)] = 0; 

pr^'^zz + IJ''&zz'i'{r + R cos 1^) — h\i^) + prRi)'lsmi^ — mgrsui{(p + d) = 0. 

(22) 

Here we have simplified the writing by defining the parameter 

p := Ks- mv^ ; (23) 

this will have a relevant role in the following. 

In the following it will be of interest to consider the Lagrangian Ctw produc- 
ing the travelling wave equations (22) as associated Euler-Lagrange equations; 
this is also obtained by restricting the Lagrangian (13) to the space of functions 
satisfying (20). It results 

= (1/2) [mr2/3(^)v2 + MR^v"^ - Kt 

-Kypi^)] dl + (1/2) [r^mv^ - K,)] 

+ [((m + M)gRcos'd + mgrcos{ip + 1?)) + h{^p)] . 



5 Series expansion for TW solutions 

A way to obtain the simple pendulums chain (and the sine-Gordon equation in 
the continuum limit) from our model is to let the length of second pendulums go 
to zero. Note that in this case our set of parameters becomes redundant, as the 
masses m and M coincide in space, so that only the total mass M = m + M \s 
relevant; and similarly for the coupling constants Kg and Kt, with total coupling 
strength K = Ks + Kt. 

If the double pendulums chain is seen as a (singular) perturbation of the 
simple pendulums one, one is naturally led to look for travelling wave solutions 
as perturbations of the standard sine-Gordon solitons. 

We will thus now look for solutions to the (22) equations in the form of a 
series expansion in a small parameter e, 

i) = ^o+£^i+e'^'&2 + .... , ^25) 

ip = ipo + epi+ £^(^2 + 

We will correspondingly also expand in the same parameter the geometrical pa- 
rameters and the masses appearing in our model, and also allow for modification 
of the speed by expanding it as well: 

r = To + sri + e'^r2 + ... , 

R — A — ro — eri — e^r2 — ... ; 

m = mo + enii + 6^1712 + ... , (26) 
M = M — mo — smi — £^m2 — ... ; 

V = Vo+£Vi+ £^W2 + ••• ; 
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note that A and M represent the total length and total mass of each double 
pendulum in the chain, which are kept constant as e is varied. In other words, 
we are now rescaling the model with s. 

We would like this is done so that for e ^ the double pendulums chain 
reduces to a simple pendulums chain, with same total length A and total mass 
M; this requires ro = 0. The choice of mo is at this point inessential (in the 
limit we have a, single pendulum with mass M = M + m), but for the sake of 
simplicity we will choose mo = 0. 

It should also be noted that once the double pendulum reduces to a simple 
one, the presence of two coupling constants Kt and Kg is actually redundant, 
as they intervene (in the limit, i.e. for £ — > 0) exactly in the same interaction: 
only K := Kg + Kf is relevant. Thus we will also write 



Kt = eki + e^k2 + ... , 
Ks = K - eki- £^^2 



(27) 



Finally, we note that for e ^ and hence for the single pendulums chain, 
the angle ip makes no sense; in order to avoid any paradoxical behavior, we 
would require that it gets frozen to zero for e ^ 1, i.e. 

-^0 = 0. (28) 

Summarizing our discussion, the series expansion we adopt are as follows: 



1? 


= 'do + s{59) 


= i^o + e-di +£^^2 + •••• , 




= e(5^) 


= eipi +s^(p2 + ■■■■ ; 


r 


= e{Sr) 


= £ri + £^r2 + ... , 


R 


= A- e{5r) 


= A - eri - e^r2 - ■■■ ; 


m 


= e{S'm) 


= £mi + £^m2 + ... , 


M 


= M-e{dm) 


= M — mg — £mi — £^m2 — 


V 


= 1^0 + £{Sv) 


= f + £-^1 + £^^2 + ••• ; 


Kt 


= e{5k) 


= £fci + £^^2 + ... , 


Ks 


= K-s{6K) 


= K - ski - £^^2 - ... . 



(29) 



We will then insert the series expansions (29) in the Euler-Lagrange equa- 
tions for TW solutions (22). It will be convenient to define 



{Mv^ - K) . (30) 



5.1 Terms of order zero 



At order zero, the second equation in (22) is identically satisfied (it just reduces 
to /i'(0) — 0, which always holds for h even), while the first reads 

i?o" = sin(t?o) := -K sin(i?o) • (31) 
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This reproduces as obvious by construction - the sine-Gordon equations 
for the single pendulums chain. Its solution is 

^}o{z) = 4 arctan [exp(K 2:)] . (32) 

Derivatives of this function are readily computed, and we have 

^0 = 2^^, V = -2^lli^^. (33) 

COSh(K2;) COSh^(K22;) 

These will appear in the higher order term of the expansion of (22); in these we 
also find terms like sin('i?o) and cos('!?o)j which in view of (32) - and with some 
simple algebra - are given by 

. , „ sinhfKx) , smh^iKx) — 1 
sm(t?o) = -2 J. \ , cos(t?o) = \/ , • (34) 

cosh [KX) cosh [KX) 

5.2 Terms of order one 

At order e, and using the equation at order £°, the first of (22) reads 

^1" = [^2^]-i \^(^{A^ -l)ki+2AMvo{Avi-nvo))^o"+ 
+ {AM g cos ^o)'&i - Mgn sm-&o ■ 



(35) 



Using the solution for ^0 and its consequences recalled above, this simplifies to 
X \a^K {cosh{2Kz) - 3) T?! + 4 B sinh(K2;) 



(36) 



where we have written 

B := (^A^Kri-Mg{{l-A^)ki+2AMvo{riVo-Avi))j . (37) 

As for the order s terms in the second of (22), these yield 

It should be stressed that these two equations, (36) and (38), are qualita- 
tively different: in facts, while the first is a (second order) differential equation 
for -di, which depends on -do, the second is an algebraic identity which directly 
determines yi. 
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5.3 Higher order terms 

Proceeding with the analysis at higher and higher orders, we would of course 
obtain more and more involved explicit equations. However, the main feature 
displayed at order one will be present at higher orders fc > 1 as well. 
That is 'dk will be determined by a differential equation of the form 

^fe" = Ffe[^?o, •••,i?fe-i; vii •••><^fc-i] + Gfe[i?o,-,'!?fe-i;</?i,-,</3fe-i] 

(omitting for the sake of brevity the dependence on parameters). 

This can be seen as the equation of motion for a particle of unit mass, whose 
position is described by '&k{t), in the time-dependent potential 

Vk{'dk,t) := -{l/2)Fk{Wkf - Gkm'&k) , (39) 

where we have used the fact '!?j(t), <fij{t) with j < k can be determined by 
solving equations at lower orders so to write^ 

Fk{t) := Fk[Mt),-,^k-i{t);Mt),-,Vk-i{t)] . (40) 

As for the term ipk, this will not be determined by a differential equation, 
but rather by an algebraic relation, i.e. 

ipk = Hk[&o, ■■■,'&k-i;V>i, ■■■,^k-i] ■ (41) 

E.g., at order s^, and using the results from the analysis of terms of order 
and , we get for this algebraic relation* 

V2 = [9/{2A^tl''h"m]x 

X \A'^KMrisech^{Kz){cos{2Kz) - 3){h"{0)fi9i 

+ AKMri[ki - A^ki + A{jlri + 2Muo(''i^^o - Av{))]x 

xsech(Kz)tanh(Kz)(/i"(0))2 (42) 

+ A^sm{§Q)[-2pi''{Q){-AK^MTi{^{f 

+ (fciMn - KMr2 + minCil + Mv^))h"{0)) 

- gK^M^rlsm{'&o)h"'{0)] . 

It should be noted that the algebraic - rather than differential - character 
of the relation (41) is due, in algebraic terms, to the choice ro = (which 
also means we are considering a singular perturbation [16]). This makes that 
the terms with derivatives of ipk will not appear in the 0{e'') terms of the 
expansion, while ipk itself does; hence we have an algebraic relation and not a 
differential equation. 



^This obviously requires to solve the equations order by order, as always in perturbative 
approaches. 

*The only interest of such an involved formula is to show that the algebraic relation can 
be explicitly determined at each order, by standard recursive computations. 



11 



In physical terms, the ehoiee tq = corresponds indeed to considering the 
double pendulums chain as a perturbation of the single pendulums one, the 
perturbation parameter being related to the length of second pendulums. 

In facts, let us substitute according to the condensate form of (29) into (24): 
we write the resulting expansion as 

£t„ = £o + e£i + £^£2 + ^(e^) . (43) 

We will not set (po = 0. At first orders we have the expressions reported in 
Appendix B 

The complete explicit form of these is not so relevant; the important fact is 
that - as can be checked explicitly by the formulas in Appendix B - we have 

dCk/d^'k = ; (44) 

hence the corresponding Euler-Lagrange equation (no sum on k here and below) 

dCk d dCk 



difk dt dip'^ 



= (45) 



reduces to the identity 



Ak := ^ = . (46) 

These equations show that as anticipated in the Introduction to any 
order of the perturbative expansion the Lagrangian £fe depends on (pk but is 
independent of the momentum conjugated to (pk- that is, the (pk can be con- 
sidered as an auxiliary field, entering in the Lagrangian only algebraically (and 
not difi'erentially). We express this fact by saying that the field </? is a pertur- 
batively auxiliary field. We stress that ip is not an auxiliary field in the full 
Lagrangian, i.c;. its auxiliary c;liaracter arises as a consequence of the (singular) 
perturbative expansion we considered. 

In facts, the above explicit expressions for £o,£i,£2 provide the expression 
for the first few identities Ak, which all reduce to h'{ipo) = 0. 

If we look at the Euler-Lagrange equations 

._ 9£k d dCk _ , . 

t-k,j ■= - -t:-^, — = U . (47) 

dipk-j dt dip'^_j 

these reduce to the algebraic identities seen above. 
In particular, we have 

5i,o = AKri{^'^cosipo + 2{i9'ofsinipo)+ipih"{ipo)=0: 
£2,1 = AKni-d'^ cos ^n + 2{-d'o)^ sin ^o)+Vih"{'Po)=0; 
£2,0 = l-Krl^o + A(i^(2(^iri(i?^)2 + r^^'^ + 
-n{ki + miVo)'dQ)] cos (po 
-[2Krl{^'of + A{2kini^'of - 2Kr2{^'of + Kr^i^^ - 4i^rii?^,t?'i 

+2miriUo(??o)^)] sin ipQ 
+n{Krii2^'^ + ^'^) - gmi sin((^o + ^0)) + (l/2)¥>?/i"'((^o) 
+h"{'po)'P2 = 0. 
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6 Discussion and conclusions 



We have considered a model corresponding to a chain of double pendulums 
(with a constraining potential limiting excursions of second pendulums), with 
first-neighbor harmonic interactions and possibly subject to an external poten- 
tial. This can, in an appropriate limit, be considered as a perturbation of the 
usual simple pc^ndulums chain with harmonic interactions, possibly subject to 
an external potential, leading to (discrete) sine-Gordon solitons. 

We dealt with these models in the continuum limit: in this case the field 
equations we obtain describe two interacting fields (possibly in an external po- 
tential) a topological field Q{x,t) and a non-topological one ^{x,t). 

Looking for travelling wave solutions - which include in particular simple 
solitons - with velocity v, we obtain a reduction to two coupled second order 
ODEs for z?(z) and (p{z), where z = {x — vt). 

These equations posses a Hamiltonian structure and as such they can be 
studied using conservation of energy; this reduces the effective dynamics - for 
given initial data - to a three dimensional manifold. However, the lack of a 
second constant of motion prevents integrability and hence obtaining a general 
solution for the dynamics. 

We argued then that our model can be seen, as mentioned above, as a per- 
turbation of the standard simple pendulums model leading in the continuum 
approximation to sine-Gordon equations. The perturbation parameter e should 
lead not only field expansion, but also appropriate expansions for the geomet- 
rical and dynamical parameters of the model. We would thus expect that ~ for 
small e - the solitons of our model can be expressed as perturbations of the 
familiar sine-Gordon solitons. 

By performing explicitly the perturbative expansion and solving the resulting 
equation at first order, we have shown that this is the case. However, in this 
procedure we obtained a rather unexpected feature: that is, the non-topological 
field is determined by the topologic;al one via an algebraic relation ~- not a 
differential equation - and thus has the role of a slaved field, while the topological 
field takes the role of master field. 

This feature is due to the appearance of a perturbatively auxiliary field; that 
is, when we expand perturbatively our model around the "do soliton, the field 
ip = eifi +e'^(p2 + ••• is an auxiliary field, determined by algebraic equations, up 
to any desired order - but it is not an auxiliary field for the full dynamics. 

We have already noted (see appendix A for details) that the simple pendulum 
system can be also obtained, acting on the confining potential, as a dynamical 
limit of the double pendulums chain. Thus, the double pendulums chain can 
be reduced to a single pendulums chain in two conceptually (and physically) 
independent ways: the dynamical way discussed in Ref. [11] and the singular 
limit discussed in this paper. Although leading to the same simple physical 
system, the two different reductions have peculiar features that make them 
rather different. Therefore they may be used in different contexts to model the 
dynamics of realistic systems such as molecular chains. 

The most striking difference between the two simple pendulums reductions 
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is related to the breaking of boost (Lorcntz) symmetry. Owing to the presence 
of two different sound speeds (related to the presence of two couphng constants 
Kt,Ks), the lagrangian (13) is not invariant under Lorentz transformation [11]. 
However, performing the e ^ limit, at the zeroth order in the perturbation 
theory we recover Lorentz invariance and the speed of the sine- Gordon soliton 
is not fixed. The recovering of the boost symmetry is essentially a consequence 
of the fact that in the limit e ^ only the coupling K = Kg + Ki and the 
mass M = m + M (not the single couplings or single masses) are relevant. It 

follows that in this limit we have just one (and not two) speed of sound K/M. 
Obviously, the boost symmetry is broken by higher orders in the perturbative e- 
expansion, but perturbative solution can be still found for any value of the speed 
of the travelling wave. Conversely, in the case of the dynamical single pendulum 
reduction the Lorentz symmetry remains broken also after the reduction and the 
soliton speed is fixed (see appendix A). Moreover, in this later case it is very 
difficult to find perturbative solutions of the double pendulums chain and one 
has to resort to numerical calculations [1]. 

In conclusion, the dynamical reduction to a single pendulums chain provide 
us with a nice mechanism to fix the speed of sine-Gordon solitons, can be applied 
in presence of a strong confining potential but is not suitable for a perturbative 
evaluation of the solutions of the double pendulums model. On the other hand, 
the singular e-expansion discussed in this paper provide us with a nice perturba- 
tive framework for evaluating the general solutions of the model, can be applied 
in presence of weak confining potentials but does not allow for a soliton speed 
fixing mechanism. 

Appendix A. 

The dynamical simple pendulums limit 

Now we note that if we force ^{x,t) = - and hence (^(z) = - in our model, 
then we are actually considering a chain of simple pendulums; in the continuum 
limit, this leads to a sine-Gordon equation. 

The constraint '^{x,t) = can be accommodated, in our setting, by acting 
on the confining potential Uc, i.e. on h: this should be made stronger and 
stronger hence we speak of a dynamical limit, as opposed to the geometrical 
one considered above - and the maximum angle (pQ smaller and smaller. 

In the limit^ — > 0+ and h"{Q) — > +oo, we expect to recover the solitons 
of the sine-Gordon equation. Note that in this limit the coupling constants Kg 
and Kf actually refer to the same interaction; we will thus take Kt = 0, also in 
order to make comparison with the sine-Gordon case more immediate. 

If we force i^{z) = 0, and use h'{0) = 0, Kt = 0, then the equations (22) 

^Wc stress that this is a singular limit, as the hmlting model (0o = 0) has a smaller number 
of degrees of freedom than any other case {cpo ^ 0). 
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reduce to 

[fi{2rR + r^ + R'^{l+p))-KsR'^p]'&^;,+gm{r + R + Rp)sm^ = 0; 
pr{r + R)t)zz + gmr sini? = . 

(Al) 

Each of these two equations provide a different determination for i?^^; the 
compatibiUty condition for them is obtained requiring the two determinations 
coincide. When r 7^ (if r = the second equation is identically satisfied) the 
compatibility condition reads 

M = -KsR/r . {A.2) 

This provides a unique determination for p and hence, in view of (23), for i>^; 
more precisely, this yields 



V mr 

Thus, the com,patibility condition selects uniquely the speed of the travelling 
wave obtained in the dynamical simple pendulums limit. 

This phenomenon is not peculiar to our model, but common to a large class 
of multi-fields soliton models, as discussed elsewhere [11]. 

Appendix B. 

Expansion of the travelling wave Lagrangian 

In this Appendix wc provide explicit formulas for the first terms in the pcrtur- 
bative expansion (43) of the travelling wave Lagrangian (24). These are: 

Co = {A/2) \A{Mvi - K){'d'o)^ - 2gMcos'do - 2h{0) 

Ci = [(1/2)(A2 - l)ki + An{K - Mvl) + A'^MvoVi - AKn cos(po i'&o)'^ 

- A\K - Mvi){'d'o&^) - AKn cos^o(i?oVo) 

- 5iM(a'!?isin^?o -ncos'i?o)^o - Wi'Po)]^! ; 

£2 = {Ay2){Mvi - K){d'^f - {l/2)Krl{^'of " AKn cos^oW'i) 
+ (1/2)[(A2 - l)k2 -2n{Aki - Kn) + 2AKr2 + Mrlv^ - 2AMr2vl 

-4.AMnvoVi + A^Mvl + 2A^MvaV2 

+2{Akin + Krl - 2AKr2 + 2Aminvl) cosLpo + 2AKnVi sin (yJo](^o)^ 
+Mgr2 cosi?o + (l/2)AgMi9f cos-i?o - gmm cos((/3o + ^0) 
-gM{rii9i - A1&2) sin^o - ¥'2/i'(<^o) - (l/2)(p?/i"(^o) _ 

+ [A'^{Mvi - K)d'^ + ((^2 - l)fci + 2AKn - 2AMnvl + 2A^MvoVi 
~2AKn cos ipo)-d'-^ + {{AKm + Krl - AKr2 + Ammvl) cos ipo 

+AKn<pi sm.^pQ - Krl)\'&'Q . 

Those expressions allow to check explicitly (at first orders) that (44) holds, 
and hence (45) reduce to (46). 
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